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Abstract: Standard deviation measurements of intensity profiles of 
stationary single fluorescent molecules are useful for studying axial local- 
ization, molecular orientation, and a fluorescence imaging system's spatial 
resolution. Here we report on the analysis of the precision of standard 
deviation measurements of intensity profiles of single fluorescent molecules 
imaged using an EMCCD camera. We have developed an analytical expres- 
sion for the standard deviation measurement error of a single image which 
is a function of the total number of detected photons, the background photon 
noise, and the camera pixel size. The theoretical results agree well with 
the experimental, simulation, and numerical integration results. Using this 
expression, we show that single-molecule standard deviation measurements 
offer nanometer precision for a large range of experimental parameters. 
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1. Introduction 

Single-molecule-fluorescence imaging has been a powerful tool in particle localization and 
tracking studies [1-4]. In single-molecule imaging, the fluorescence intensity profile of a point 
light source is called a point spread function (PSF). While the PSF is described by an Airy func- 
tion, it is, in practice, approximated by a Gaussian. A Gaussian fit to the PSF of a stationary 
single fluorophore has two fitting parameters: centroid and standard deviation (SD). The cen- 
troid is the center of the PSF, and the SD is its width. While the centroid determines the lateral 
position of the particle in the imaging plane, the standard deviation determines its axial position 
[5-7] and orientation [8-10], as well as the spatial resolution of the fluorescence imaging system 
[11]. 

Error analysis provides the precision for a physical measurement, and is essential for valida- 
tion of the method used. While error analysis of single-molecule PSF centroid measurements 
has provided the precision for lateral localization measurements [1], which enabled differenti- 
ation of various biological mechanisms (such as the walking mechanisms of myosin V on actin 
[2]), PSF SD measurement error analysis will provide the precision in the following applica- 
tions: (1) single-molecule axial position measurements, where the SD of a single molecule's 
PSF increases with the defocusing distance [5-7]; (2) single-molecule orientation measure- 
ments, where at different molecular orientations with respect to the imaging plane, the molecule 
exhibits an elliptical-shaped Gaussian PSF with a SD that changes in both lateral directions de- 
pending on its orientation [8-10]; and (3) characterization of a single-molecule-fluorescence 
imaging system, where the measured SD of an imaged single fluorophore determines whether 
the imaging system is diffraction-limited [11]. 

In contrast to the precision of centroid measurements which has been extensively investigated 
and applied to many systems [1,2, 12], the precision of SD measurements of single-molecule 
PSFs has not been evaluated. Here we report SD measurement error studies of immobile-single- 
molecule PSFs using analytical calculation, numerical integration, simulation, and experimen- 
tal measurements. As with centroid analysis, the precision of SD measurements is affected by 
the experimental settings of a finite number of photons per PSF N, the standard deviation of 



the background noise Ob, and the camera's finite pixel size a. We have obtained an analyti- 
cal expression for the PSF SD measurement error as a function of these parameters. Our SD 
measurements have achieved nanometer resolution for a wide range of experimental condi- 
tions. This expression for the SD measurement error will provide confidence in determining a 
particle's axial position and molecular orientation from measurements using a single-molecule 
imaging system of known resolution. 

2. Theory 

2.1. Formulating SD measurement error, As, by % 2 minimization 

The term "standard deviation" was introduced by Pearson in his 1894 mathematical study of 
evolution [13] and characterized further in the following years [14]. For different collections 
of size N of randomly selected data from a common distribution with a theoretical standard 
deviation sq, the error associated with the SD measurement of each collection is ^/sq/2N, 
as first calculated by Pearl in 1908 [15]. The same expression was derived more recently by 
Taylor using a different method [16]. In this article, we derive the SD error for a PSF, which is 
a collection of photons from a common distribution emitted by a point light source. We include 
the additional experimental effects of photon count fluctuation per PSF, background noise, and 
camera pixelation in our study. 

We utilized the method developed by Bobroff [17] and subsequently used for centroid error 
analysis by Thompson, Larson, and Webb [1] to derive the error associated with SD measure- 
ments of single fluorophores. The approach uses Chi-square statistics to estimate the error as- 
sociated with fitting of experimental data to expected theoretical values. In order to maintain 
consistency in notation for single-molecule tracking studies, we will retain many of the same 
notations used in Ref. [1]. Below, we derive the analytical solution to the PSF SD error as a 
function of N, a, and Ob (b in prior studies) beginning with one dimension and extending to two 
dimensions. 

In ID least squares fitting of the intensity profile of an immobile single fluorophore, % 2 {s) 
is proportional to the sum of squared errors between the observed photon count at pixel i, y„ 
and the expected photon count Nj(x,s), of a PSF. Here x and s are the measured position and 
SD of the PSF, respectively, while xq and sq are the true location and the theoretical SD of the 
molecule: 

z 2 w=E^, a) 

i "i, photon 

where O^ pnoton is the expected photon count uncertainty at pixel i without accounting for 
photon-to-camera count conversion (described in the following section). In this article, we 
emphasize the SD error and assume that the location measurement errors are negligible, i.e. 
x = xq (Appendix B shows that the codependence of localization and SD errors vanishes). For 
simplicity, Nj(xo,s) is denoted as Nj in this article unless otherwise specified. 

There are two sources for Oi^hoton at pixel i: one is the Poisson-distributed photon shot noise 
of the PSF where the variance is the mean expected photon count of the pixel, A 7 ,-, and the other 
is the SD of the background noise, Ob, expressed in photons. The variances of the two sources 
add to yield 

tfphoton=Ni + °Z- (2) 

The deviation of s from so, As =s — sq, is obtained by setting dx 2 (s)/ds to 0, expanding Nj 



about sq, and keeping the first order term in As: 
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where AT- is the derivative of A 7 ,- with respect to 5 evaluated at so, and Ay,- = iVj(jcO)*o) — By 
squaring Eq. (4) we obtain the mean squared value of As, 

Li(N/ 2 / a f photon) 

The root mean square of As, As rms , is the PSF SD error that we calculate in this article. Appendix 
A shows the detailed derivation of Eq. (5) from d% 2 {s) jds = 0. 

2.2. Modifying Oi iP h ton to include camera count conversion effects 

When an EMCCD (Electron Multiplying Charge Coupled Device) camera is used in imaging 
single fluorescent molecules, the detected pixel reading is in camera counts. In converting from 
camera counts to photon counts, an additional variance in O^ pnoton appears. Below we derive 
the uncertainty in photon counts, a,, to use in place of Oi^photon m Eq. (5) for experiments where 
EMCCD camera count conversions are involved. 

An EMCCD camera amplifies the detected photons by converting each photon to a distribu- 
tion of photoelectrons through many multiplication stages. At the final stage, one photon yields 
a distribution of camera counts (equivalent to the last stage photoelectron counts) with a dis- 
tribution function /(«*) [18], 

/(«*) = I exp(-nVM), (6) 

where n* is the camera counts in the distribution and M is the photon multiplication factor of 
the camera. Here we use * to denote camera counts in order to differentiate from photon counts. 
The n* distribution has a mean of M and a variance of M 2 . 

At pixel i, the PSF photon count distribution is described by a Poisson distribution with the 
variance being equal to the mean. Each photon at the pixel contributes two terms to the pixel's 
camera count variance: the mean photon shot noise variance M 2 (variance of a single photon 
count, which is one, multiplied by the square of the multiplication factor), and the photon- 
to-camera count conversion variance M 2 . The total camera count variance contributed by one 
photon is 2M 2 ; therefore, a mean of Ni photons yields a camera count variance of 2NjM 2 . This 
variance agrees with the expression in Ref. [19] where the variance in camera counts a 2 ut camera , 
is related to the variance in photon counts (7 (n p noton , by an excess noise factor F , 
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for EMCCD cameras with a large number of multiplication stages. 

Fluorescence from buffer, diffusing molecules in the solution, and camera counts from elec- 
tronic readout and thermal noise constitutes the total background photon count at pixel i, with 



a variance of of and a mean of (b) . The total background variance in camera counts is the sum 
of the background count variance <J%M 2 , and the variance introduced by the average number of 
background photons, (b), each with a variance of M 2 : (o£ + {b))M 2 . 

Summing the PSF and the background contributions, the total camera count variance at pixel 
;' is 

a* 2 = 2NiM 2 + (a 2 h + (b) )M 2 . (8) 
When expressed in photon counts, 

a 2 = af/M 2 = 2N l + a 2 + (b). (9) 

Revising Eq. (5) with the modified a, we have 

2.3. Expressing As in photon counts 

To evaluate Eq. (10) in ID, we use a normalized Gaussian distribution 

JV,-=^Uxp(-H 2 /2. 2 ), (11) 
\J2its 



where we set the location of the PSF to be at xq — for simplicity and without loss in generality. 
We approximate the pixel summation in Eq. (10) by an integral going from negative to positive 
infinity, and we estimate ((As) 2 ) at the two extrema of of: the high photon count regime where 
a? + (b) can be neglected, and the high background noise regime where 2Ni can be neglected. 
In the high photon count regime, 

<(As) 2 > = |, (12) 
and in the high background noise regime, 

((As)2)= vWM 03) 

An alternative derivation of Eq. (12) is presented in Ref. [16], although the photon-to-camera 
count conversion variance was not included and thus ((As) 2 ) — s 2 J2N. The total ID ({As) 2 ) is 
the sum of Eqs. (12) and (13) (without the pixelation effect discussed below) 

«^> = f + V 1iS +W ' 

The method of approximating ((As) 2 ) by summing these results for both extrema of of is 
validated by numerical calculation results shown in Fig. 2, and is in accordance with Ref. [1]. 
We now calculate the effect of camera pixelation on ((As) 2 ). Each photon in a PSF is associated 
with two variances with respect to the centroid. One is the mean variance of the PSF, s§, and 
the other is due to the fact that each photon is further binned into a pixel that has an intensity 
profile described by a uniform distribution with a width corresponding to the pixel size a. The 
variance of this distribution is a 2 / 12. Thus, the total variance of a photon due to pixelation is 
the sum of the two, 

4 2 =*o+^- (15) 
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Under experimental conditions, the measured s should be (s§ + a 2 /12) and for theoretical 
formulations, the expected SD of a PSF should include the pixelation effect. We have verified 
that sj 2 increases with a according to Eq. (15) by simulation. Plugging Eq. (15) into Eq. (14) 
we have for ID 

, (16) 

u ' 1 N 3aN 2 

Extending the ID ({As) 2 ) calculation to 2D where s x #, which for the remainder of this article, 
represents the SD in either the x or y direction of the imaging plane, and sq x and SQ y are the 
theoretical SD values in the x and y directions, respectively, 

„ Af *v _ 4 + i , 16^(4+f^) 3/2 (j+f^) 1/2 K 2 + W) m , 
U* 1 *))- N + 3a 2 N 2 ■ 

The derivation of Eq. (17) is provided in Appendix C. 

A more accurate estimation of ((As x ^ y ) 2 ) can be obtained by numerically integrating Eq. (10), 
incorporating the transition region between the high photon count and the high background 
noise regimes. The numerical integration results are shown in Fig. 2 to be consistently higher 
than the analytical calculation results by w 15%. 

3. Methods 

3.1. Experimental setup 

Single-molecule imaging was performed using a Nikon Eclipse TE2000-S inverted microscope 
(Nikon, Melville, NY) attached to an iXon back-illuminated EMCCD camera (DV897ECS-B V, 
Andor Technology, Belfast, Northern Ireland). Prism-type Total Internal Reflection Fluores- 
cence (TIRF) microscopy was used to excite the fluorophores with a linearly polarized 532 nm 
laser line (I70C-SPECTRUM Argon/Krypton laser, Coherent Inc., Santa Clara, CA) focused to 
a 40 jUm x 20 jUm region on fused-silica surfaces (Hoya Corporation USA, San Jose, CA). The 
incident angle at the fused-silica water interface was 64° with respect to the normal. The laser 
was pulsed with illumination intervals between 1 ms and 500 ms and excitation intensity be- 
tween 0.3 kW/cm 2 and 2.6 kW/cm 2 . By combining laser power and pulsing interval variations 
we obtained 50 to 3000 photons per PSF. A Nikon 100X TIRF objective (Nikon, 1.45 NA, oil 
immersion) was used in combination with a 2X expansion lens, giving a pixel size of 79 nm. 

At focus, the PSF image generated by a point light source with a mean emission wavelength 
of 580 nm and symmetric polarization has a full width at half-maximum (FWHM) of m A/2NA 
= 580 nm/2.9 w 200 nm and theoretical s = FWHM/2.35 w 85 nm. Including the pixelation 
effect [Eq. (15)], the measured PSF SD Sq x0 ^, for our imaging system should be 88 nm. Due 
to random fluctuations in the emission polarization direction of streptavidin-Cy3 molecules 
attached to surfaces [8] and variations in focus between each measurement, we observed a 
range of Sq x values from 90 nm to 140 nm. 

Single streptavidin-Cy3 molecules (SA1010, Invitrogen, Carlsbad, CA; 530/10 excitation, 
580/60 emission) were immobilized on fused-silica surfaces by depositing 6 /J.I of 0.04 nM 
streptavidin-Cy3 powder dissolved in 0.5X TBE buffer (45 mM Tris, 45 mM Boric Acid, 1 mM 
EDTA, pH 7.0). A coverslip flattened the droplet and its edges were sealed with nail polish. The 
fused-silica chips were cleaned using oxygen plasma before use. We inspected for possible sur- 
face fluorescence contaminations by imaging the TBE buffer alone; no impurities were found 
on either the fused-silica surface or in the buffer. The immobilization of the adsorbed molecules 
was verified by centroid vs time measurements. 



3.2. Data acquisition and selection 

Typical movies were obtained by synchronizing the onset of camera exposure with laser illumi- 
nation for different intervals. The gain levels of the camera were adjusted such that none of the 
pixels of a PSF reached the saturation level of the camera. For the initial step, streptavidin-Cy3 
monomers were first selected in ImageJ (NTH, Bethesda, MD) by examining the fluorescence 
time traces of the molecules for a single bleaching step [20]. For a selected monomer, the in- 
tensity values for 25 x 25 pixels centered at the molecule were recorded. The center 15x15 
pixels of the PSF were used for 2D Gaussian fitting with peripheral pixels used for background 
analysis. 

The intensity values of the selected molecules were first converted to photon counts (see the 
following section) and then fitted to the following 2D Gaussian function using a least squares 
curve fitting algorithm (lsqcurvefit) provided by MATLAB (The Mathworks, Natick, MA): 

/(^)=/oex P (-^-^) + W, (18) 

where fo was the amplitude and (b) was the mean background value. A background pixel's total 
count is the sum of the floor, electronic readout noise, and background fluorescence counts. For 
the (b) in this article, the floor value, determined by the lowest background pixel value, has 
already been subtracted. With this fitting, the PSF's SD values in both the x and y directions, its 
measured location (xo,yoX an d the image's mean background value were obtained. 

The selected streptavidin-Cy3 monomers were further characterized to satisfy the following 
conditions used for SD error analysis. (1) No stage drift detected by using centroid vs time 
measurements. Stage drift introduces additional blur to each single-molecule PSF and thus 
affects the measured SD values. (2) A minimum of 75 valid PSF images, each with a photon 
count N that fluctuated less than 20% from the experimental mean (N), of the monomer. The 
PSF N count restriction is necessary for precise SD error analysis at N by using a statistically 

sufficient number of PSFs with consistent N. (3) PSFs with signal-to-noise ratios (Iq / \J h + of) 

larger than 2.5, where Iq is the peak PSF photon count (total photon count minus (b)) and of 
is the background variance in photons. (4) Mean (s x ) and (s y ) obtained by Gaussian fitting of 
the s x and s y distributions of all valid images did not differ by more than 10 nm, or ±5% of the 
mean SD value to minimize polarization effects of Cy3. (5) The mean SD values (s XJ ) were 
between 95 nm and 135 nm to minimize defocusing effects. These constraints on s x and s y are 
necessary for obtaining the expression for As rms , as a function of N, with minimal variations in 
the other parameters. 

3.3. Photon gain calibration 

To convert from a pixel's camera count to photons, the camera count value was divided by M. 
In order to obtain M for each experimental setting, the center nine pixel values of the PSF were 
evaluated if the molecule's average signal-to-noise ratio was greater than 3. Conversely, when 
the signal-to-noise ratio was less than 3, only the center pixel was used for calculation as the 
adjacent pixels contained too few photons for statistically accurate calculations. According to 
Eq.(9), 

M = {af-af)/im) -(b*)), (19) 

where (N*) and a* are the Gaussian fitted mean and standard deviation of the measured camera 
count distribution of pixel i, respectively. Here (N*) is the mean camera count that includes 
background fluorescence and electronic noise counts. For each image, the average M for all 
nine center pixels of the PSF (or just the center pixel for PSFs with low signal-to-noise ratios) 



was obtained. Then, the average M values for all fitted PSFs in a movie were again averaged to 
obtain the camera's multiplication factor for this movie. 

In order to verify that our method of calculating M is correct, we have simulated PSFs with 
low and high N, including background photon noise and confirmed that Ni follows a Pois- 
son distribution which approaches a Gaussian at high N with variance Nj + a?. Including the 
photon-to-camera count conversion variance [Eq. (6)] in the simulation, we verified Eq. (19). 

3.4. PSF and background simulations 

Single-fluorescent-molecule PSFs were generated using the Gaussian random number genera- 
tor in MATLAB. For Fig. 2, the SQx.Oy of each simulated PSF was determined by the experi- 
mental means (s x _ y ) . The observed fluctuation in the number of photons N, was incorporated. 
The generated photons of each PSF were binned into 15x15 pixels with a pixel size of 79 nm. 
Then each photon count in a pixel was converted into camera count using Eq. 6 with a M value 
of one. Random background photons at each pixel were generated using the corresponding 
experimental background distribution function. Although the exact experimental background 
distributions were used for the simulations, the numerical integrations and analytical calcula- 
tions were computed using the theoretical variance and the mean of all background counts, a% 
and (b) , respectively, rather than their fitted values. The background counts are primarily drawn 
from two types of distributions: a full Gaussian with a high mean or a truncated Gaussian with 
a low mean (Fig. 1C), depending on the background fluorescence level of each specific experi- 
ment. The final simulated PSFs with background noise were fitted to a 2D Gaussian [Eq. (18)] 
to obtain the centroid and SD values of the PSF. 

For each simulated As x ^ rms data point, 1000 iterations (2000 iterations for Fig. 3) were per- 
formed and the Gaussian fitted SDs of the s x>y distributions were the simulated As x ^ rms results. 

4. Results 

We report our study of 2D As x>rms using four different methods: (1) experimental measurements, 
(2) simulations, (3) numerical integrations of Eq. (10), and (4) analytical calculations using Eq. 
(17). 

Figure 1 A shows a set of single streptavidin-Cy3 molecule images with an increasing number 
of detected photons N. These molecules have similar mean SD (s x ) values of 1 10 nm, 111 nm, 
and 107 nm, respectively. In order to demonstrate the decreasing SD error with increasing N, 
each representative image was chosen such that the 2D SD value was the sum of the mean 
SD (s x ), and one standard deviation of the molecule's s x distribution As X:rms (SD,„ Mi , e = (s x ) + 
As*,™*)- To clearly illustrate the change in the SD error, which is measured as the PSF SD 
minus (s x ), the ID intensity profiles of the PSFs are plotted in Fig. IB as opposed to their 2D 
intensity profiles for clarity. The ID intensity values were obtained by averaging transverse 
pixel intensity values of the PSF at each longitudinal pixel i. It is evident that the widths of the 
ID Gaussian fits decrease with increasing N. The measured 2D SD imc!S , e values deviate from 
their respective means, (s x ) values, by 10.3 nm, 7.2 nm, and 2.7 nm. Again as expected, when 
N increases, the 2D SD error decreases. 

Figure 1C presents the background distributions associated with the molecules in Fig. 1A. 
The background distribution function resembles either a full or a truncated Gaussian, depending 
on experimental settings. The corresponding oj, and (b) values for these images are 0.92 and 
2.75 photons, 0.81 and 2.07 photons, and 1.48 and 1.97 photons, respectively. The lines are fits 
to the distribution and the histograms represent our simulated results. 

Figure 2 shows As x rms obtained by using experimental measurements, simulations, numeri- 
cal integrations, and analytical calculations. Each experimental As x jrms data point is the standard 
deviation of the s x distribution for a single streptavidin-Cy3 monomer. A simulation was per- 
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Fig. 1 . (A) Representative images with increasing A' of 151, 393, and 1891 photons of single 
streptavidin-Cy3 molecules. It is evident that the "blurriness" of the molecules decreases 
with increasing N. The "blurriness" is defined to be the measured SD of the image and is 
the sum of the mean SD of the molecule's PSF and its SD error. (B) ID intensity profiles 
(circles) of the molecules in (A) and their Gaussian fits (lines). The respective ID SD values 
are 195.4 nm, 140.5 nm, and 1 10.9 nm, and the respective deviations of the 2D SD values of 
the images from their means are 10.3 nm, 7.2 nm, and 2.7 nm. As expected, these deviations 
from the mean decrease with increasing N. The scale bar is 500 nm. (C) Background count 
distributions (circles) for the three molecules in (A) and their fits (lines). The histograms 
are simulated background distributions which reproduce those observed experimentally. 
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Fig. 2. Comparing As x ,rms vs N obtained by using four different methods: experimental 
measurements (solid squares), simulations (circles), numerical integrations (crosses), and 
analytical calculations (dashed line). Each experimental As xrms data point is the SD from 
the Gaussian fit to the s x distribution of a single streptavidin-Cy3 monomer. For each data 
point, its experimental N and background distributions were used for simulation, and its 
experimental (N), (s x , Y ), o"/,, and (b) values were used for the numerical integrations and 
analytical calculations. The experimental data are on average 57% higher than the analyti- 
cal calculation data. 



formed for each experimental data point. The parameters were based upon experimental results 
including fluctuations in a PSF's total detected photons, background distribution, and the so^Oy 
values determined by the mean experimental (s x ,y) after subtracting for the pixelation effect 
[Eq. (15)]. The finite bandwidth of the emission filter was also taken into consideration by sim- 
ulating each photon as being drawn from a PSF whose width is varied according to a Gaussian 
distribution centered about so x fl y (with SD of 2 nm). Numerical integrations and analytical cal- 
culations used the same (N), so x fly, C7&, and (b) as those in the corresponding experimental data 
point. For all N, the numerically integrated As x ,,-„„- results are pa 15% higher than the theoretical 
results while the experimental results are sa 57% higher, and the simulations agree well with 
the experimental results. 

The above results are for our pixel size of 79 nm. For different experimental settings the 
pixel size will vary and affect As xrms . Figure 3 shows As xrms vs a/ so studied by simulations 
and analytical calculations using sq x = soy = sq — 120 nm, N = 500 photons, Oj, = 1 photon, 
and (b) — 4 photons. The generated photons of each PSF were binned into 19x19 pixels and 
subsequently converted into camera counts following the same procedure described above for 
Fig. 2. As a /so increases, there is an initial decline in As xslns until rising at a /so ~ 0.73. Beyond 
a/ so 0.73, As x>m s increases slightly and then continues the decline again at a /so ~ 1.18. 
This decline after a/so ~ 1.18 disagrees with theory, which suggests an increase in As xrms 
beyond the theoretical minimum of (a/ so) 4 — — at a /so — 1.18 (vertical dashed line). 

47T( CT 2 + (i.)) + 

The overall decreasing As x rms trend after the theoretical minimum occurs because when the 
pixel size increases, the measured PSF SD is increasingly affected by the width of the pixel 
and approaches the SD of the pixel; thus, variations among measured SD values decrease. 
Eventually, at sufficiently large pixel sizes where the whole PSF is contained within one pixel, 




Fig. 3. Adeems vs a/ so studied by simulations (circles) and analytical calculations (dashed 
line; theoretical results shifted up by 57%) for N = 500 photons, .so = 120 nm, Cj, = I 
photon, and (b) = 4 photons. In these simulations, there were no fluctuations in N and 
s 0x = *Qy = s 0- Note that the general trend illustrates that As x ,rms decreases with increasing 
pixel size. For a/ so < 1.18, the simulated As xrms results agree with the theoretical results 
multiplied by (1.23 + 0.42a/io) (solid line). The shifted theory line (up by 57%) crosses 
the simulation line at a/so ~ 0.80 (a ~ 96 nm). Considering the random discrepancy of 
up to ±15% between simulation and experimental results in Fig. 2, the cross point lies 
within an acceptable range given our experimental pixel size of 79 nm. At a /so > 1.18, the 
experimental As x , rms results continue to decrease influenced by the increasing dominance 
of the pixel's SD. The vertical dashed line at a/so — 1-18 is where the theoretical As x>rms 
minimum occurs, determined by differentiating Eq. (17) with respect to a. 



the measured SD will be the SD of the pixel, inferred by the top-hat distribution function, 
and the measured SD error will be zero. The analytical calculation does not take this large 
pixelation effect into consideration; consequently, these results and those of the simulations 
begin to rapidly diverge. 

The simulated local As xrms minimum occurs at a/so = 0.73, rather than at the theoretical 
minimum of a /so = 1.18 due to the pixel size effect described above. Our experimental settings 
of a = 19 nm and so = 120 nm yield a /so = 0.66 and is close to the simulated As xjms minimum. 
The dashed line is the theoretical results shifted up by 57% and the solid line is the theoretical 
results multiplied by (1.23 + 0.42a/*o), yielding an excellent fit to the simulation results for 
a/ so < 1.18. We have also performed additional simulations using different parameter sets 
where the theoretical minimum always preceded the continued decline in As xrms . According to 
Fig. 3 and our other simulations, a good a/ so range for future studies should be between w 0.5 
and 1, as is usually the case. Future As Xi y >rms studies using different pixel sizes should take this 
discrepancy into account. 

Note that the simulated As Xtrms minimum at a/ So = 0.73 is different from the theoretical 
Ax nm minimum at a/so = 0.88 described in Ref. [1], and our theoretical Ax rms minimum at 
a/so = 1.10 calculated from Eq. 20, using our set of parameter values. Future studies should 
take this difference into consideration by selecting an optimal pixel size. 



5. Discussion and Extensions 



Here we discuss four issues: (1) causes for discrepancies between results obtained using dif- 
ferent methods; (2) modifications to the centroid measurement error developed by Thompson, 
Larson, and Webb [1] to include the EMCCD camera photon conversion effects; (3) relation 
between SD error and the error of the measured quantities associated with each of the afore- 
mentioned applications; and (4) methods to determine the SD error As xjsms , for dimeric fiuo- 
rophores and mobile molecules in future studies. 

5.1. Causes for discrepancies 

Numerical integration results are consistently higher than the analytical results by 15%, while 
simulation results are higher than analytical results by 57% for all N. There are a number of 
reasons for these discrepancies: (1) The analytical As x ^ rms result [Eq. (17)] is obtained by eval- 
uating Eq. (10) for the two limiting cases of of at the high photon count and high background 
noise regimes. The intermediate regime is absent and thus the numerical integration and simula- 
tion results are larger. (2) When Ni is expanded about sq, the higher order terms were neglected 
[Eq. (4)]. (3) In the As rms calculation (Appendix A), the Ni distribution function is assumed to 
be a Gaussian for all pixels of the PSF [Eq. (23)]. This assumption will only be statistically ac- 
curate for center pixels of PSFs with high N. For peripheral pixels, especially for PSFs with low 
N, the Ni distribution function approaches a Poisson with a low mean, rather than a Gaussian. 
These different Nj distributions, which have been verified by simulation, were not considered 
in the analytical calculations. (4) In simulations, we attempted to model the background count 
distribution exactly, whereas in numerical integrations and analytical calculations, the shape 
of the background count distribution was not considered, and therefore did not influence the 
results. 

In summary, the analytical calculation of the SD measurement error expressed in Eq. (17) is 
a reasonable approximation for a large range of experimental parameters. When the 57% dif- 
ference is corrected for, the expression is in excellent agreement with our experimental results. 
Future studies using this formula should be aware of the limitations and be sure to include this 
57% difference from underestimation of the true error for similar a/so values. 

5.2. Modifications to centroid error analysis 

The PSF centroid error expression developed by Thompson, Larson, and Webb [1] did not take 
the photon-to-camera count conversion variance into consideration. Additionally, the theoreti- 
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cal standard deviation so, should be modified to include the pixelation effect («Q + a 2 /12) , 
with respect to both directions. We have modified the PSF centroid measurement error to be 

, rAx rt_ 2(4+g) , 8^(.L+f^) 3/2 (4+n) 1/2 K 2 + W) f2m 

U*X))- N + ^ . (ZU) 

This theoretical expression for the centroid measurement error underestimates the experimental 
results by 42%. 

5.3. Interpreting As X y, rms in SD measurement applications 

With regards to the three applications of SD measurements discussed in this article, the SD 
measurement error of a single image can be translated into the precisions associated with each 
of the application's measured quantities. For the direct translation, SD error is the uncertainty of 
an imaging system's measured resolution; for the indirect translations, the precisions for axial 
localization and molecular orientation measurements can be expressed by SD error. 



For characterization of an imaging system's spatial resolution, if the system is diffraction- 
limited, the SD of the imaging system should fall between the measured SD of the PSF from a 
visible point light source ± the SD error [11,21]. Thus, SD measurement error directly provides 
the precision for quantifying an imaging system's resolution. 

For axial localization studies, it has been previously shown that the SD of the PSF of a 
molecule located at a distance z away from the focal plane can be expressed as [5-7] 

( z 2 \ l/2 

s(z)=so(l + jp) , (21) 

where D s=y 400 nm is the imaging depth of a typical single-molecule imaging system. Conse- 
quently, by error propagation, the precision in the SD measurement of a single image, As rms , 
can be used to determine the localization error associated with the molecule's axial position, 
Az: 

Az = — tjt Aw (z) . (22) 

For molecular orientation studies, the polarized PSF, for a range of orientations, has an ellip- 
tical intensity profile that can be fit by a 2D Gaussian with different standard deviations in the 
x and y directions [8-10]. When an expression relating s x and s y to the orientation is developed, 
the error in measuring s x and s y , once again by error propagation, can be used to calculate an 
error associated with the reported orientation of the molecule. 

5.4. As x . y,rms calculation for future SD measurement applications 

In addition to stationary single molecules, SD measurements can be used in future applications 
to investigate molecules such as stationary dimers or moving fluorophores. We are currently 
exploring these two areas of interest: (1) two sub-diffraction limit, separated molecules labeled 
with identical fluorophores that exhibits a combined PSF with a SD that increases with their 
separation [22]; (2) a moving molecule (i.e. a freely-diffusing fluorophore) which produces a 
blurred image given a finite exposure time, whereby the measured SD of the resulting intensity 
profile can be used to study the dynamic properties of the molecule, particularly its diffusion 
coefficient [23]. With modification, the method for estimating the SD error of stationary single 
molecules in this article can be extended to these two cases. For these studies, the Ni distribution 
function at each pixel may be different from the Gaussian assumption for stationary molecules 
in Eq. (1 1). A new Ni distribution function for each specific case can be obtained and a new of 
formula [Eq. (9)] can be derived. Using the new Ni distribution function and of, the SD error 
for these cases can be obtained following the same procedure outlined in the theory (Sec. 2). 

6. Conclusion 

In this article we report the precision analysis for SD measurements of single-fluorescent- 
molecule intensity profiles. Our analytical expression of the PSF SD error allows for proper 
quantification of the precision associated with determination of the imaging system's resolution 
and both axial localization and molecular orientation measurements of single molecules. Fur- 
thermore, we propose additional studies to characterize multiple fluorophores and examine the 
diffusive properties of mobile molecules by evaluating the measured SDs of their correspond- 
ing intensity profile to known precision. When our theoretical framework is extended to these 
studies, SD analysis will be advanced into a powerful tool for single-molecule-fluorescence 
imaging studies. 



Appendix 



Here we present the complete derivation of Eq. (5). We first obtain a probability distribution 
function for y,-. At large N of a few hundred photons, the y, probability distribution function 
at each of the center nine pixels of the PSF is a Gaussian, while at the peripheral pixels, the 
y,- probability distribution function is better approximated by a Poisson with a low mean. Here 
we assume that our N is significantly larger than 100 photons and the y,- probability distribution 
functions for all PSF pixels are Gaussian functions 

1 __/ Ay? 



* = 73ii exp l-#J' (23) 

where Ay,- = N t (x ,so)- yi and of is rfphoton as in Eq. (1). For Gaussian distributed y,-, we have 

(Ay,) = 0, (24a) 
((Ay,) 2 ) - of. (24b) 

Starting from Eq. (1) and taking a derivative with respect to s, 

d% 2 (s) _ v d (yj-N,) 2 _ v 2(y l -N l )(y l -N i y^-(y l -N l ) 2 -2a l g' i 

Setting the above equation to zero, we find 

V 2(y l -N l )(y,-N,y „ (y, - N,) 2 ■ 2g,g/ 

I ~ f =1 ~ t ■ <*> 

We can simplify Eq. (26) using the following terms: 



y t -N,{s) = y t - (Ni(s ) + Nj As) = - Ay i-Nl As, (27a) 

{Ji-Ni)' = —N[, (27b) 

of = 2Ni{s)+2ol = 2(N i (so)+N' i As)+2oj l , (27c) 

lOio'i = 2N[. {21 A) 

Inserting Eqs. (27a)-(27d) into Eq. (26), we obtain 

^ -2(Ay,+jV»Hv;) = y (Ayj+NjAsf-lNl 

i of i °i 

^ ^ (Ay 2 +2Ay,jV;A.)-2jV; 

i 

Moving As to the left-hand side, 



This equation is Eq. (3). Neglecting the Ay,/ a, 2 term, we get Eq. (4). 



We now take the mean square of Eq. 4. Note that the average is meant to apply to y, only, so 
we have 
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((As) 2 ) = "' " = • . (30) 

For two different pixels, their distributions are independent, so (AyiAyj) = 5y((Ay ( ) 2 ) = of 
[see Eq. (24b)]. This gives us Eq. (5). 

B. 

Here we calculate the codependence of Ax and As in ID and show that it vanishes. Thus the 
assumption that x (or s) is fixed when taking a partial derivative of Ni with respect to s (or x) is 
valid. 

Expanding Ni(x,s) about (x ,s ) to first order => A/j = A^xq^o) + A*7srU + A^^|. So . 

Setting = = and solving for As and Ax, the error associated with the s measurement 
is 
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The cross product term ( 3£i^7ifo j = 0> an d arriving at Eq. (5), ((As) 2 ) 
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1.-57; 
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In Appendix C we calculate the 2D ((As x 2 )). In 2D, the expected counts at pixel i, j is given by 

^^^(jg.-<*£\, ,32, 

where we assume that the PSF is centered at zero. Taking the derivative of Ni with respect to s x 
and evaluating at sq x , 

((As 2 )> = ^^. (33) 



it 

1 ds x 



Next, we approximate the summation by an integral where ;' and j are continuous from negative 
to positive infinity. There are two limits to the approximation, one being the high photon count 

limit and the other being the high background noise limit. At the high photon count limit, 

2 

((As*) 2 ) = 'jf- after taking the photon-to-camera count conversion variance into consideration. 
At the high background noise limit, ((A^) 2 ) = lto Woy(g|,+W) ^ Adding ^ e two terms together 

1 /2 

and replacing SQ X fi y by (sq x 0y + a 2 /12^ to incorporate the pixelation effect, we arrive at Eq. 
(17). 
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